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Abstract 

We consider a spaceborne gravitational wave (GW) detector formed by three spacecraft in a 
linear array, coherently exchanging laser beams and using the data combinations of time-delay 
interferometry (TDI) . We previously showed how the measured time series of Doppler shifts in the 
six one-way laser links between spacecraft pairs in a general unequal-arm triangular configuration 
can be combined, using TDI, to exactly cancel the otherwise overwhelming phase noise of the 
lasers while retaining sensitivity to GWs. Here we apply TDI, unfolding the general triangular 
configuration, to the special case of a linear array of three spacecraft. We show that such an 
array ("SyZyGy") has, compared with an equilateral triangle GW detector of the same scale, 
degraded (but non-zero) sensitivity at low-frequencies (f << c/ (array size)) but similar peak and 
high-frequency sensitivities to GWs. We develop the GW and noise responses, noting that in this 
geometrical case only one TDI combination is GW-sensitive showing a relatively simple 6-pulse 
response to an incident GW burst. Sensitivity curves are presented for SyZyGys having various 
arm-lengths. A number of technical simplifications result from the linear configuration. These 
include: only one faceted (e.g. cubical) proof mass per spacecraft, intra-spacecraft laser metrology 
needed only at the central spacecraft, placement in an appropriate orbit can reduce Doppler drifts 
so that no laser beam modulation is required for ultra-stable oscillator noise calibration, and 
little or no time-dependent articulation of the telescopes to maintain pointing. Because SyZyGy's 
sensitivity falls off more sharply at low-frequency than that of an equilateral triangular array, it 
may be most useful for GW observations in the band between those of ground-based interferometers 
(~ 10 - 2000 Hz) and LISA (~ 10~ 4 - 10" 1 Hz). SyZyGy with ~ 1 light-second scale could, for the 
same instrumental assumptions as LISA, make observations in this intermediate frequency GW 
band with 5a sensitivity to sinusoidal waves ~ 2.5 x 10~ 23 in a year's integration. 
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I. INTRODUCTION 



When coherent optical or microwave wave trains are interchanged between three or more 
freely orbiting spacecraft, and phase or frequency differences of all incoming and outgoing 
beams recorded, there is a closure that allows subsequent elimination of all source fluc- 
tuations ("laser noises") and so achieves sensitivity to propagating gravitational waves at 
much lower levels set by secondary noises. In doing this, the precise propagation times 
between sources and readout locations enter the data reduction algorithms, but equality of 
arm lengths is not necessary. We have called this scheme "time delay interferometry" or 
TDI, Q and subsequently analyzed several configurations and data types needed for further 
elimination of the most serious secondary phase noise sources, viz. those due to spacecraft 



non-geodesic motion ( "optical bench noise" ) and especially to fluctuations of the u 



tra-stable 



local oscillators, or frequency standards, used in fringe tracking ("USO noise") . Fringe 
tracking is necessary when the lasers in the configuration have widely differing frequencies, 
and/or when the configuration does not move rigidly so that emitted and received signals 
are Doppler shifted during slow spacecraft separation changes. 

The LISA space gravitational wave antenna will use TDI, and is to be roughly equilateral, 
which maximizes the overall response to low frequency gravitational waves over a wide band. 
Reduction of secondary noise, and simplicity of system architecture, nevertheless suggest 
consideration also of other schemes; here we consider the most extreme possible application 
of TDI, viz. to a linear orbiting array, or, in astronomical terms, three spacecraft flying in 

syzygy 1 - 

We have previously considered a linear array, or parallel beam interferometric detector 
of gravitational waves, in work showing that with multiple bounces a split antiparallel laser 
beam will amplify a gravitational wave signal while laser fluctuations are cancelled conven- 
tionally by precise control of equal arm lengths [6] . We now see the multiple readout scheme 
of TDI as much preferable, and equal spacing of the three spacecraft is no longer required. 
The equations for combining the heterodyne data streams in TDI can be taken directly from 



1 An amusing coincidence is that Dhurandhar, Nayak and Vinet |(| have pointed out that underlying the 
data combinations of TDI is the mathematics of polynomial rings in three variables, and that the structure 
we found is called by algebraists a "module of syzygies" ! SyZyGy for this reason again would seem to be 
an apt acronym for the concept of this paper. 
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our previous papers, which had general triangular configuration; in this paper they will be 
used in the linear limit. 

The (nearly) linear formation of SyZyGy has some advantages and disadvantages over 
the nominal equilateral triangular formation planned for the LISA mission. The primary 
disadvantages of the linear formation are less sensitivity at the low frequencies (for a given 
spacecraft separation distance) and sensitivity to only one linear combination of GW po- 
larizations. The advantages arise from the potential for a more stable formation geometry. 
The nominal LISA triangular formation involves spacecraft in orbits with semi-major axis 
a, inclination I, and eccentricity e related by e = J/3 1 / 2 and I = L/(2a) 3, Ol- This choice 
of orbital elements produces a rotating triangular configuration with arm lengths that are 
constant to ~ e 2 L. The time variation of the arm lengths causes a variable Doppler shift on 
the laser beam signals, a variable direction from each spacecraft to the other two spacecraft, 
and a variable angle between the spacecraft velocity and the direction to the other space- 
craft. The (almost) linear array might be formed with three spacecraft moving in a single 
circular orbit about the sun. In the ideal case (ignoring perturbations from other bodies) 
the distance between spacecraft and the directions from each spacecraft to the others will 
remain constant. The constant distance between spacecraft removes the Doppler shift on 
the laser signals. The constant directions means that no mechanism is needed to account for 
angle variation. The constant angle also means that a single faceted (e.g. cubical) test mass 
can be used in each spacecraft, instead of the two in the nominal LISA configuration. (In the 
case of spherical test masses, the constant angle between the velocity direction and direction 
to other spacecraft removes some surface irregularity errors that would affect a triangular 
formation with spherical test masses.) The constant angle between the velocity direction and 
direction to other spacecraft means that the point-ahead and point-back angles are constant. 
(These angles arise from the need to point lasers to hit the location the other spacecraft will 
be at after a one-way light time.) The linear combination needs less spacecraft propulsion 
to reach the desired orbits (for a given location relative to Earth and arm length L) since 
no plane change is needed. No laser metrology is needed for the proof masses relative to the 
optical benches at the outlying spacecraft (e.g., data we have denoted Zi 2 , z 32 , at spacecraft 
2). (At the center spacecraft (spacecraft 1) the intra-spacecraft metrology data denoted z 2 \ 
(or — %l) in the TDI equations must still be taken to eliminate phase fluctuations due to 
vibrations of the optical benches.) 
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The procedure for elimination of USO noise introduces considerable complication in the 
LISA system, inasmuch as additional calibrating data must be taken and exchanged if present 
generation flight-qualified USO's are used. The USO noise enters the TDI combinations 
multiplied by the internally generated fringe tracking frequencies u>ij that are required when 
there are inter-spacecraft Doppler drifts and/or offsets of the laser center frequencies. Con- 
trol of relative spacecraft motions will automatically be achieved by spacing SyZyGy along a 
single nearly-circular solar orbit, and laser frequency offsets can be limited either by atomic 
line stabilization or by use of optical transponders [jj. Thus no wide-band fringe tracking 
data, USO calibration, and associated modulation of the laser beams are needed by SyZyGy, 
as will be shown in the next section. 

In the linear case the gravitational wave antenna patterns of all the TDI data combination 
become simple and identical (so-called 6-pulse response). We give this pattern explicitly in 
Section 2. The pattern will of course sweep over the sky as the spacecraft move along 
their common orbit. We apply the TDI data combinations in the linear geometry, and 
show how the only data combinations that are sensitive to GWs have the same signal 
response as the combination X, the unequal arm Michelson interferometer. We introduce 
two new independent combinations that have no gravitational wave response, denoted p and 
cr; they can be used to monitor system noises. All previously identified laser-noise-free data 
combinations are given in terms of X, p and a. 

In Section 3 we derive the SyZyGy response functions for the remaining secondary noises: 
optical bench vibration, proof mass residual accelerations, optical path noises at the pho- 
todetector readouts, and USO noise entering in to the heterodyne measurements. In Section 
4 we give numerical calculations and plots for these. Averaging over the sky and polarization 
states, we finally plot the resulting gravitational wave sensitivities to incident sine waves, us- 
ing one year observation time and several SyZyGy scale sizes. Compared with an equilateral 
LISA, a SyZyGy of the same scale has similar best and high-frequency gravitational wave 
sensitivities. However the sensitivity is poorer at low-frequencies, suggesting applicability 
in the frequency band between those of ground-based interferometers and LISA. 
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II. ANALYSIS 

We use the notation and results of 

as 

. A straight interferometric array, shown in Figure 
1, is obtained by opening up the general triangular configuration of Figure 2 of reference j^, 
so that spacecraft 1 is in the center and L\ = L 2 + L 3 . We denote the linear orientation by 
the unit vector h — hi — —fi2 = —n^. On each of the spacecraft the usual pair of optical 
benches can be fused into one, so we put V* — t^. At the outlying spacecraft 2 and 3 both 
incoming beams will be referenced to the same face of a single proof mass. At the center 
spacecraft, number 1, TDI will require use of two opposed faces of a single proof mass. For 
all we will have v* = V{ and the optical bench and proof mass motion components that enter 
are Vi ■ h and v-i ■ h, respectively. The intra-spacecraft data z^ at the three proof masses 
consequently satisfy 232 = -212,-213 = -223, and z 2 \ = —Z31. The first four of these then drop 
out of all the laser-noise-free data combinations and so in fact need not be recorded. The 
last two remain, and so one of them needs to be taken, which amounts to implementing 
precision laser metrology solely on the central spacecraft. 

When all these specializations are substituted into the data combinations derived for gen- 
eral configurations, it is found that they degenerate into essentially one, six pulse, response 
to incident gravitational waves (we use X) plus two others that respond only to system 
noises (we call them p and cr). All three will further be USO-noise-free to acceptable limits 
if the laser frequency offsets are limited by Doppler effects due to relatively stable spacecraft 
motion in the SyZyGy linear orbital configuration. 

Specifically, the new laser-noise-free combinations we use for a linear array are 

X — 2/32;3'22' ~ 2/23;233' + 2/31;22' _ 2/21;33' + 2/23;2 — 2/32;3' + 2/21 — 2/31 — z 21 ~ z 21;ll' + ^21;22' + ^21;33' 
P = 2/13 - 2/23 - 2/31;2' + ^31;2' 
& = 2/12 - 2/32 - 2/21;3 + ^21;3 

where a semicolon subscripted index i is understood to mean not only time delay by Li but 
also multiplication by the Doppler factor (1 — Li/c), which in this case can be dropped 2 . 
Note how interchanging indices 2 and 3 takes p and a into one another. 

2 In more general work 0, Q the TDI combinations have been modified for the case of a moving array where 
time delays such as Li from spacecraft 2 to spacecraft 3 differ from L 1 from spacecraft 3 to spacecraft 2, 
etc. One then has for SyZyGy Li = L 2 + L 3 and L 1 = L 2 + £3. The primes on the time-delay subscripts 
in these equations for X, p, and a have been inserted for that case, but we do not insist on that generality 
in any of the subsequent discussion. 
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The previously used TDI basic combinations in this linear array limit (and without mo- 
tion) are given in terms of X, p, and a: 
a = X + p. 2 - cr ;3 

(3 = -<T + P;1 

7 = P ~ <?;1 
C = P;3 - <T ;2 

Proof mass and shot noise response functions for X and p (or a) are given in Section 4. 

III. NOISE TRANSFER FUNCTIONS 

The proof-mass-plus-optical-bench assembly and laser beam paths for the central space- 
craft is show schematically in Figure 2. The photodetectors that generate the time series 
l/2i) J/3i, and 2 2 i are as indicated. The outgoing light from spacecraft 1 to spacecraft 2 is 
routed from laser 1 on the optical bench using mirrors and beam splitters; this beam does 
not interact with spacecraft l's proof mass. Conversely, the incoming light beam from space- 
craft 2 is first bounced off the nearest face of proof mass 1 before being reflected onto the 
photodetector, where it is mixed with light from laser 1. This time series is y^i in Figure 
2. Analogously, the time series 1/21 is producted from light transmitted to and received from 
spacecraft 3. The relative motion of the optical bench and proof mass is monitored by mea- 
suring the Doppler shift of the light coming directly from the laser with the light bounced 
off one face of the proof mass (this forms the metrology time series z 21 = —231, as indicated 
in Figure 2. 

Figure 3 shows the light paths on spacecraft 3. The outgoing light beam from spacecraft 
3 to spacecraft 1 and 2 is routed from laser 3 on spacecraft 3's optical bench using mirrors 
and beamsplitters; again, this beam does not interact with spacecraft 3's proof mass. The 
two incoming light beams from spacecraft 1 and 2 are bounced off the proof mass before 
bing reflected onto photodetectors where they are mixed with light from laser 3. The time 
series produced are y 23 and yi 3 , as shown. The configuration at spacecraft 2 is analogous, 
producing y 12 and y 32 . 2/13 and y i2 enter a, p, and a, but not X. 

The equations for the noises entering the and z 2 \ = —Z31 can be developed from 
Figures 2 and 3. Consider first the situation on spacecraft 1. The photodetector producing 
y 3 i (moving with velocity Vi) reads that time series by mixing the beam originating from 
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the optical bench on spacecraft 2 sent in direction n 3 = —n (laser noise C 2 and optical bench 
motion V 2 , delayed by propagation along L 3 ), after one bounce off the proof mass (vi), with 
the local laser light (Ci). The time series y 2 \ is produced analogously The z 2 ± measurement 
is from light originating at the laser (Ci, Vj), bounced off the left face of the proof mass (vi), 
and mixed with the direct laser light (Ci). For spacecraft 2 and 3, y 32 , and 1/23 are produced 
similarly. Metrology time series on the end spacecraft are not required. 

We will see below that there is no gravitational wave contribution to the data p and a. 
It may however be desirable to measure these to monitor system responses and secondary 
noises. In that case we would also need the data y± 2 and yx3- 

The yij and the required metrology time series z 2 i (= — Z31), including gravitational wave 
signals and shot noises, can be developed from the diagrams and by consulting the Appendix 
of Tinto et al. M : 



^0 2/2i = v 3 (l + ^3,2 - n 2 ■ V 3>2 )A 2 (1 + 2n 2 ■ v x - h 2 ■ Vi) 
+ Ci) - W2i(l + Qi) + z^o 2/fr + u o V° 2 \ Patk 



(1) 



vq 2/31 = ^2(1 + C 2 , 3 ' + n 3 • V 2i 3')A3(l - 2n 3 • 0i + n 3 • Vi) 
-i/i(l + C x ) - w 3 i(l + Qi) + i/ 2/fr + ^0 2/31*' Pa<?i 



(2) 



1/0^21 = 2^n 3 • (vi - Vi) = -i/ 23i (3) 

where A$ = (1 — L»/c). Other are determined from equations (1) and (2) by cyclic index 
permutation. 

Inserting these into the equations for a and p we can verify that the laser phase/frequency 
noises cancel and that the optical bench noises and USO noises are, for reasonable parameters 
of SyZyGy, suppressed to acceptable levels. Noting that A 2 A 3 ~ Ai and u 23 — u X3 = 
z/iA 2 — 1/2 A 1 ~ —^31, and its permutations, we then obtain the following responses to 
proof-mass, USO, and optical path noises in X and p or a. 
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j^noise j^proof mass _|_ j^USO _|_ j^opt. path 

= -2(z/i/V )(n • V! - n ■ v^u) + 2(z/ 2 /z/ )(n ■ v 2 - 3 - n ■ ^322) 

+2(1/3/1/0) (n • V 3t2 - fl ■ V 3 -233) 

-2(uj 31 /u )(n ■ vi -h-Vi-22) -2(uj 21 /u )(n-v 1 - n • v 1;33 ) 
+2(u 32 /vo)(n ■ v 2 - 3 - n ■ v 2;322 ) + 2(cj 23 /z/ ) (n • v 3 - 2 - n ■ v 3 . 233 ) 
+(^31/^0) (<5i ~ Qi&i) ~ (^21/^0) (Qi ~ Qijsa) 

+ (^32/^o) (<52;3 - <?2;12) ~ (Wus/^oXQ^ - <?3;13) 

. opt. path opt. path . opt. path opt. path 
+l/32;322 — 2/23;233 + ^31;22 — 2/21;33 

■ opt. path opt. path . opt. path opt. path /a \ 

+i/23;2 ~~ #32:3 + S/21 ~ Vsi l 4 J 



pnoise pproof mass _|_ U SO _|_ p°Vi- path 

= 2(u 31 / v Q ) [n-v 3 -n- v 1;2 ] 

-(^3l/t / o)(<53 - Qlfl) 

. opt. path opt. path \ opt. path /r\ 
+2/31 ~ V23 ~ ZA 2^31;2 V 3 ) 

Under the assumptions that all noises are independent, that optical path noise only depends 
on the distance between spacecraft pairs, and that for proof mass noise in X terms with 
coefficients ujij/iy are negligible compared with terms having coefficients of order unity, the 
spectra of the noises in X and p are: 



O ! t\ rtproof mass . nUSO . nopt. path 

ox{j) = b x + b x + b x 

~ 16[sin 2 (27r/L 1 ) + sin 2 (27r/L 2 ) + sin 2 (27r/L 3 )] 1 S proo/ mass 

+ [4(^ 2 2 + 0J 2 31 )/vl sin 2 (27r/L 2 ) + 4(w 2 3 + c^V^o sin 2 (27r fL 3 ) 

-8(^ 3 i)/z/ 2 cos(27r/(L 2 -L 3 )) sin(27r/L 2 ) sm(2irfL 3 )]S uso 

+8 sin 2 (27r/L 3 )5 2 f • path + 8 sm 2 (2irfL 2 )S 3 f ' path (6) 
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<j c<proof mass _|_ cjUSO _|_ £topt. path 

_ of. ,2 /,,2\ oproof mass , r>/, ,2 /,,2\ qUSO 

, C*°pt- po>th i <^°pt- pot/i i gopt. po&h 

13 23 31 v / 



IV. GRAVITATIONAL WAVE SIGNAL RESPONSE 

The gravitational wave signal response in TDI combination X was given for a general 
triangle in equation (24) of reference The SyZyGy configuration leads to considerable 
simplification. Taking the limit as the general triangle is distorted to a line, using L\ = 
L 2 + L 3 , and remembering that h — fi\ — —h 2 = —n 3 , we obtain a six pulse response with 
a single amplitude depending only on the orientation of the linear array with respect to the 
source direction: 



X 9W = 2 k ■ n [$(t - L 2 + L 2 k ■ h) - 

+ ^>(t - L 3 - L 3 k ■ n) - if?(t - L x - L 2 - L 3 k ■ h) 

+ #(i -Lx-Li- L 3 ) - #(i - L x - L 3 + L 2 k ■ h)} (8) 

where * is the Rational waveband * is (see QQ) 

U = 2 l-(fc-n) 2 ( } 

for a transverse traceless gravitational wave h(t) = [h + (t) e + + h x (t) e x ]. Note that X gw 
is zero for waves incident at right angles to SyZyGy (k • n = 0) and for waves parallel or 
antiparallel with it (k ■ h = +1 or —1). 

Although SyZyGy will not operate exclusively in the long-wavelength limit (LWL - 
2 r nfL\ << 1), analytical results at low-frequency are useful. Expanding the gravitational 
wave h in series of spatial derivatives and substituting into X gw gives 



X 9W -> -iiL x L 2 L 3 h ■ h'" ■ n (10) 

where ji = k ■ h. 
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V. GRAVITATIONAL WAVE SENSITIVITY 



Figure 4 shows root-mean-square X gw , averaged over polarization and the celestial sphere 
jl|, ^, for 3 armlength cases: (a) L 2 = L 3 = 16 light-seconds, (b) L 2 = 7, L 3 = 9 light- 
seconds, and (c) L 2 = 0.7, L 3 = 0.9 light-seconds. Notice that in the LWL X gw oc f 3 , as 
expected from the LWL expansion above. 

Figure 5 shows the noise power spectra of X, using proof mass and optical path noises 
only, for the same armlength cases as Figure 4. (USO can be made negligible compared to 
the other noises through trajectory control: if the maximum relative speed is kept below 
0.2 m/sec-either through trajectory design or by periodically using thrusters to adjust the 
maximum Doppler shift-USO noise will not contribute.) In constructing the curves in Figure 

5 we used the same instrumental spectra as in our previous sensitivity analysis of LISA. The 
proof masses were assumed to have independent noises with one-sided acceleration spectral 
density 3 x 10~ 15 m sec -2 Hz -1 / 2 (this corresponds 

SB 

to a one-sided spectral density of 
relative Doppler S pro °^ mass = 2.5 x 10~ 48 [f/l Hz]~ 2 Hz _1 . Following previous practice 
we approximately account for all optical path noise, including beam pointing noise (see 
and reference |l3(, Table 4.1), on a given link with S° ptical path = 1.8 x 10" 37 [f/1 Hz] 2 (r/5 x 
10 9 m) 2 Hz _1 , where r is the distance between spacecraft on that link. The aggregate noise 
spectrum plotted uses these spectra and the X transfer function (see equations 4 and 6). 
The GW sensitivity is proportional to the ratio of the rms noise to the rms signal. Figure 

6 shows the sensitivity of SyZyGy expressed in the conventional way: sensitivity to isotropic 
sinusoidal gravitational radiation (SNR = 5 in one-year integration), 5[Sx(f) -B] 1 ^ 2 /rms 
X 9W (f), where B = 1 cycle/year, for the three armlength cases of Figures 4 and 5. 



VI. DISCUSSION 



SyZyGy offers some technical simplifications compared with other configurations. Each 
spacecraft utilizes a single proof mass, and laser metrology is needed only at the central 
spacecraft. Because of the low relative velocities afforded by the SyZyGy orbit, the com- 
plexity needed for USO noise elimination (laser beam modulation and additional calibration 

n 

data) |5j is not required by SyZyGy. The angles between the arms are much more constant 
in time than for LISA, so time- dependent articulation of the telescopes is not required. 
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SyZyGy combinations p and a have exactly zero GW response and can thus be used with 
no ambiguity to assess on-orbit noise performance. Deliberate resetting of arm lengths (i.e., 
different arm lengths in different portions of a long mission) is easier in a linear configuration 
and could be used to tune SyZyGy's best response over its mission lifetime. 

For similar scale size, SyZyGy achieves both its best sensitivity and high-frequency per- 
formance comparable to those of LISA; at low-frequencies, SyZyGy has poorer sensitivity. 
This suggests applicability in the Fourier band between those of LISA and ground-based 
interferometers. Little attention has been given to the ~ 0.1 — 10 Hz band. The work that 



3- 



has been done suggests it to be relatively free of foreground (astrophysical) GW sources 
A SyZyGy with L\ ~ 1 light-second could test this with 5cr sensitivity ~ 2.5 x 10~ 23 in a 
year's integration. 
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FIG. 1. SyZyGy configuration. Unit vectors hi point between spacecraft pairs with the 
indicated orientations. The arrows on the Lj indicate the sense of the un-primed light times. 
Each spacecraft has one optical bench and one proof mass. 

FIG. 2. Signal routing and readouts on the central spacecraft, 1. 

FIG. 3. Signal routing and readouts on spacecraft 3, one of the end spacecraft. 

FIG. 4. RMS gravitational wave response of SyZyGy, TDI combination X, as a function 
of Fourier frequency for three cases: L 2 = L3 = 16 light-seconds; L 2 = 7, L3 = 9 light 
seconds; and L 2 = 0.7, L 3 = 0.9 light seconds. 

FIG. 5. One-sided noise spectra of X as a function of Fourier frequency. One-sided 
spectra of proof mass noise and optical path noise are taken to be the the same as we have 
used in previous calculation of LISA sensitivity : respectively, 3 x 10 m sec Hz ' 
and 20 x 10~ 12 x (r/5 x 10 9 m) m Hz _1//2 , where r is the distance between spacecraft. Center 
frequencies of the lasers are assumed to be equal and a maximum Doppler shift (longest 
arm case) of 0.2 m/sec is assumed. Solid lines: spectra of X including shot and proof mass 
noises. Dashed lines: contribution to X noise spectrum from an uncalibrated USO having 
one-sided spectrum of fractional frequency fluctuations 8 x 10~ 27 /[f/lHz] (see text). 

FIG. 6. Gravitational wave sensitivity of X for SyZyGy (SNR = 5 in a one year inte- 
gration, sky averaged, sinusoidal signals). Arm lengths, gravitational wave responses, and 
noise spectra are as in Figures 4-5. 
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